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URSA MAJOR: A MISSING LOW-MASS CDM HALO? 
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ABSTRACT 

The recently discovered Ursa Major dwarf spheroidal (dSph) galaxy candidate is about five to eight 
times less luminous than the faintest previously known dSphs And IX, Draco, and UMi. In this Letter, 
we present velocity measurements of seven color-magnitude selected UMajor candidate stars. Two of 
them are apparent non-members based on metallicity and velocity, and the remaining five stars yield a 
systemic heliocentric velocity of v = —52.45 ±4.27 kms" 1 and a central line of sight velocity dispersion 



of (v 2 ) 1 2 = 9.3i**2 7 kms -1 , with 95% confidence that (v 2 ) 1 '" > 6.5 kms" 1 . Assuming that UMajor is 
in dynamical equilibrium, it is clearly dark matter dominated, and cannot be a purely stellar system like 
a globular cluster. It has an inferred central mass-to-light ratio of M/L ~ 500 Mq/L© and, based on 
our studies of other dSphs, may possess a much larger total mass to light ratio. UMajor is unexpectedly 
massive for its low luminosity - indeed, UMajor appears to be the most dark-matter dominated galaxy 
yet discovered. The presence of so much dark matter in UMajor immediately suggests that it may be 
a member of the missing population of low-mass galaxies predicted by the Cold Dark Matter (CDM) 
paradigm. Given the weak correlation between dSph mass and luminosity, it is entirely likely that a 
population of dark dwarfs surrounds our Galaxy. 

Subject headings: galaxies: individual: Ursa Major dwarf - galaxies: kinematics and dynamics - Local 
Group - dark matter - celestial mechanics, stellar dynamics 
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1. INTRODUCTION 

All of the Local Group dwarf spheroidal (dSph) galax- 
ies have velocity dispersions much larger than expected 
for self-gravitating stellar systems, implying that their dy- 
namics is dominated by dark matter, with the stars being 
little more than dynamical tracers within a dark halo {e.g. 
Mateo 1998). Some dSphs have central mass to light ra- 
tios of ~ 100 M Q /L , and average mass to light ratios of 
several hundred {e.g. Kleyna et al. 2001, 2004; Wilkinson 
et al. 2004). 

Recently, Willman et al. (2005) discovered a new can- 
didate dSph in Ursa Major in a search of data from the 
Sloan Digital Sky Survey (SDSS) (Abazajian et al. 2005). 
Located about 100 kpc from the Galaxy, its half-light ra- 
dius is r\/2 ~ 250 pc, covering 7.75' on the sky. Most 
remarkably, with My ~ —6.75, it is about five times less 
luminous than the faintest previously known dSph And IX 
(Zuckcr et al. 2004; Chapman et al. 2005) and about eight 
times less luminous than the faintest Galactic dSphs UMi 
and Draco (Mateo 1998). 

In this paper, we describe precise radial velocity mea- 
surements of seven UMajor candidate stars. We compute 
UMajor's inferred central density and mass, and discuss 
how UMajor fits into the rest of the rest of the dSph pop- 
ulation within the context of standard CDM theory. 

2. DATA 

We selected candidate stars from Data Release Three 
of SDSS, picking targets within 6' of UMajor's position 
as given by Willman et al. (2005). Figure 1 shows the 
color magnitude diagram of the central 6' region, slightly 
smaller than Willman's half-light radius r 1 / 2 = 7.75'. The 
giant branch and horizontal branch are both visible, al- 



though together they contain only about 50 stars. Our 
seven targets (dark circles in Figure 1) were drawn from 
the brightest part of the giant branch, and span the mag- 
nitude range i — 17.45 to % — 18.38. Additionally, we ob- 
served the bright velocity standards HD107328, HD90861, 
and HD132737. 

We observed our stars using the upgraded HIRES (Vogt 
1994) echelle spectrograph on the Keck I telescope on the 
night of May 17, 2005. Each star was observed in one 
integration lasting 1800s. The long spectral coverage of 
the spectrograph allowed us to obtain wavelengths from 
Ha (6564A) to the red-most line of the Calcium triplet 
(8662 A). The signal to noise varied among the spectra 
from S/N = 12 to 3 per pixel, or S/N = 61 to 17 per A. 
Only the red-most Ca triplet line was near a sky line, but 
sky subtraction was generally clean even in this case. 

We extracted the spectra using the MAKEE data reduc- 
tion package for HIRES, creating flux and variance spectra 
for each echelle order. Because the latest version of MA- 
kee for the new three-chip upgrade to HIRES does not at 
present solve for the dispersion, we fit the dispersion solu- 
tion manually for the relevant echelle orders using IRAF. 
All calibration arc exposures had to be taken during the 
afternoon or morning because of persistent ghosting with 
the new chips. To verify instrument stability, we mea- 
sured the positions of sky lines in the science exposures, 
and found them to be constant within 0.2 kms -1 . 

For each of the four stellar absorption lines of interest 
(Ha and the three Ca lines), we cross-correlated a syn- 
thetic Gaussian template with the appropriate echelle or- 
der using the IRAF FXCOR package. In all cases except 
for the third Ca triplet line for the faintest UMajor star, 
we obtained a clear cross-correlation function peak. 

A potential problem with single-slit observations is that 
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mis-centering of the star on the slit will produce a veloc- 
ity offset common to all orders. To address this problem, 
we note that telluric absorption lines experience the same 
spurious velocity shift as stellar absorption lines, and can 
be used to adjust the velocity back to its correct value. 
Accordingly, we computed a velocity adjustment using tel- 
luric absorption features around 6880A. Using a template 
created from the telluric lines from one of the bright stan- 
dard stars, we used cross-correlation to compute relative 
velocity adjustments of the other stars. In all cases, the 
adjustment was less than 0.8 kms -1 . 

An additional potential source of inaccuracy is template 
mismatch; for instance, the Ha line is actually a blend of 
several different transitions, and may not be exactly at 
the template wavelength. To address this problem, we 
adjusted all of the velocities for each stellar line en masse 
relative to the second Ca line by computing the median ve- 
locity difference between the line and the second Ca line. 
This procedure introduced a shift of at most 1.6 krns" 1 . 

Next, we treated each velocity for each stellar line as an 
independent measurement, and combined them to obtain 
a final velocity and error. Finally, we shifted the entire ve- 
locity set to bring the standards into agreement with their 
published velocities, with a final scatter of 0.4 kms -1 , 
slightly larger then the intrinsic standard star uncertainty 
of 0.3 kms -1 . For each star, we obtained two error es- 
timates: o~ iraf > obtained by rescaling the nominal IRAF 
FXCOR velocity errors to give the expected x 2 ; and cr scat , 
obtained from the empirical velocity scatter among indi- 
vidual lines. For UMajor stars, the two errors are generally 
similar, but ctiraf apparently overestimates the errors for 
the bright standard stars. Table 1 shows the results for 
our seven target stars. 

In Figure 2, we show the velocities as a function of mag- 
nitude. Objects 3 and 4 are clearly outliers, separated 
about 60 kms -1 from the other five objects. 

The equivalent width of the Ca triplet can be used as a 
measure of a star's metallicity, for known surface gravity 
(dwarf vs. giant) (e.g. Armandroff & Da Costa 1991). To 
establish whether the kinematical outliers are members of 
the UMajor population, we computed the equivalent width 
of the Ca triplet lines by fitting them with a Lorcntzian 
profile plus a linear continuum, and then integrating the 
difference between the line profile fit and the continuum 
fit. We compute the uncertainties of the equivalent widths 
using a Monte-Carlo procedure: we generated simulated 
data using our best fit solution added to the empirical 
noise spectrum output by the makee reduction package, 
and fit the simulated data in the same manner as the real 
data. 

Figure 3 shows the equivalent widths for the three Ca 
lines. The kinematical outliers 3 and 4 clearly have a 
larger equivalent width for the second and third Ca lines, 
although the first line shows no difference. Object 7 pro- 
duces very poor fits because the lines are barely discernible 
above the noise, and its equivalent widths are highly sus- 
pect. Generally, only the second and third lines are used 
for metallicity determination (Armandroff & Da Costa 
1991), so the absence of an effect for the weaker first line 
is not a source of concern. However, for a fixed metallic- 
ity, the equivalent width is anti-correlated with magnitude 
with a slope about 0.6A per magnitude (Armandroff & Da 



Costa 1991), so that it appears anomalous that star 6 has 
a larger equivalent width than the brighter stars 1,2, and 
5. To be cautious, we consider the possibility that both 6 
and 7 are non-members. 

3. VELOCITY DISPERSION AND SYSTEMIC VELOCITY 

The simple RMS line of sight dispersion of the five mem- 
ber objects (1,2,5,6, and 7) is 9.1 kms -1 , to which mea- 
surement errors contribute only slightly. To compute the 
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dispersion (u 2 ) more rigorously, we assume that the ob- 
jects are drawn from a Gaussian distribution, and apply 
equation 2 of Klcyna et al. (2004), imposing an a priori 
maximum value of 50 kms -1 on the dispersion, and using 
cr SC at as the velocity uncertainty. We do not consider the 
effect of binary stars, which is generally small. Including 

all member stars, we obtain (y 2 ^ 1 ^ 2 = 9.3^ 11 2 7 kms -1 , 
where the uncertainties are lcr, defined such that 16% 
of the probability distribution is both above and below 
the error interval. If we omit object 7, then we have 

(v 2 ) 1 ' 2 = 10.4ii 6 2 9 kms -1 , and if we omit both 6 and 

7, (y 2 ^ 1 ^ 2 = 6.9^q 9 4 7 kms -1 . For these three cases, in the 
same order, we are confident with 95% certainty that the 
dispersion is greater than 6.5, 7.3, and 4.7 kms -1 . 

Assuming five genuine members with a true dispersion 
of 9.3 kms -1 , the systemic heliocentric velocity of UMajor 
is -52.5 ±4.3 kms -1 . 

4. MASS TO LIGHT RATIO 

The central mass to light ratio of a stellar system with 
an isotropic orbital distribution and constant M/L may 
be expressed as (M/L) n = 9i](v 2 } / (2-KGrhbSo), where i) 
is a luminosity distribution dependent factor always close 
to 1, rhb is the half-light radius, and So is the surface 
brightness (Richstone & Tremaine 1986). Using UMa- 
jor's approximate published total luminosity My — —6.75 
(L v ~ 4 x 1O 4 L ) and half-light radius r 1 / 2 — 250 pc 
(Willman et al. 2005), we obtain a central surface bright- 
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ness of 5*0 = 0.11 L Q pc~ 2 . Assuming that (u 2 ) = 9.3, 
one then obtains (M/L) ~ 550M Q /L Q , as one might 
expect for a system that is nearly an order of magni- 
tude fainter than known dark-matter dominated dSphs, 
but has the same dispersion. The central density is then 
po = 0.18 M Q pc -3 = 6.8GeVcm~ 3 using the above 
equation, or /?o = 0.22 Mq pc -3 using the formula /?o = 
166 (w 2 ) /r coro (Mateo 1998). Despite UMajor's extreme 
M/L, its density is within a factor of two of that of sev- 
eral other dSphs (Mateo 1998). Several important caveats 
apply, however. First, we have identified the half-light 
radius with the half-brightness radius, when in truth they 
differ. If UMajor has a Plummcr profile, for example, then 
Hib = 0.64 r 1 / 2 , and we underestimate M/L substantially. 
Next, the above formula for po/Io strictly applies only to 
systems with constant M/L, whereas a dwarf with a halo 
has varying M/L, a fact that is often ignored when com- 
puting M/L using core fitting. Thus the above values of 
(M/L) and po are intrinsically imprecise, and are useful 
primarily as relative values to compare with other dSphs 
having similarly computed qualities. 

An alternative approach to computing the mass of 
the system is the projected mass estimator (Bahcall & 
Tremaine 1981; Heisler, Tremaine, & Bahcall 1985; Evans 
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et al. 2003), whereby the mass of a system with mea- 
sured radial velocities vt at projected radii Ri is given by 
M = C/(GN) x J2iLi v iRii where C is a constant de- 
pending on the mass and light distributions of the sys- 
tem. For our five-member UMajor data, this becomes 
M = C1.5 x 10 6 M Q . Evans et al. (2003) compute C 
for the general case of measuring the mass between radii 
T'min to r max for a p oc r~ 7 tracer population in a (f> oc r~ a 
potential. Taking r m ; n = 100 pc from our innermost star 
and assuming that all the stars are within the half-light 
radius so that r max = 250 pc, and assuming that 7 = 3, 
with a halo having a flat rotation curve (a = 0), we obtain 
C w 6, again giving M « 1O 7 M and M/L ~ 500 inside 
the region observed. The mass varies from 5 x 10 6 M to 
2 x 10 7 M Q as 7 ranges from 2 to 6. 

We further note that dSphs for which radially extended 
data exists have a mean M/L within the apparent stellar 
cutoff or "tidal" radius that is nearly an order of magni- 
tude larger than the core-fit central M/L (e.g. compare 
Mateo (1998) with Kleyna et al. (2001), Wilkinson et al. 
(2004), and Kleyna et al. (2004)). If we accept the CDM 
result that halos are similar in structure, then UMajor 
might have a global M/L ~ 3000M Q /L Q . 

5. DISCUSSION 

Our stellar velocity measurements show that UMajor 
is a nearly dark, low-mass halo, almost an order of mag- 
nitude less luminous than previously known dSph galax- 
ies of similar mass. We obtain a central mass-to-light 
ratio M/L ~ 500M Q /L Q ; if this extrapolates to larger 
radii like other dSphs, UMajor may have a mean M/L ~ 
3000M Q /L Q . Hence, UMajor may be the most dark mat- 
ter dominated galaxy yet discovered. For comparison, the 
recently discovered 10 11 M & VIRGOHI 21 Hi source has 
M/L > 500M Q /L Q (Minchin et al. 2005), and the ex- 
treme dark spiral NGC 2915 has M/L w 80 (Meurer et al. 
1996). 

The CDM structure formation paradigm predicts that 
dark matter clumps into cusped halos characterized by a 
single parameter, the central density, with the smallest and 
least massive halos being the densest (Navarro, Frcnk, & 
White 1997). Problematically, CDM also predicts an or- 
der of magnitude more small Galactic satellite halos than 
are actually observed (e.g. Klypin et al. 1999), but it 
is unclear whether the theory is at fault, or whether the 
missing halos are too dark to be observed. 

High Velocity Clouds (HVCs) have been suggested as 
candidates for the missing CDM dwarf halos (Blitz et al. 
1999). This suggestion remains controversial, especially as 
the high velocity clouds may themselves be divided into 
compact and diffuse subclasses with possibly different ori- 
gins (e.g. de Heij, Braun, & Burton 2002). There is no 
clear distance determinant for the HVCs around the Milky 
Way, and, for a long time, there was no clear consensus 
as to whether the HVCs were Galactic or extra-Galactic. 
This ambiguity may have been resolved by the detection 
of a faint circumgalactic HVC population around the An- 
dromeda galaxy by Thilker et al. (2004), who argue that 
the available data are consistent with formation mecha- 
nisms via cooling flows, or with tidal debris from mergers, 
or with gaseous counterparts of the missing CDM haloes. 

Accordingly, UMajor may represent the best candidate 



for a "missing" CDM halo. Its existence raises several in- 
teresting questions. 

Why do the dSphs have a similar velocity dispersion, a 
similar central density and, presumably, a similar mass, as 
noted by Mateo et al. (1998)? It is curious that Draco, 
UMi, and UMajor, the three lowest luminosity Galactic 
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dSphs, all have the same dispersion, about (v 2 ) = 
10 kms^ 1 , and And IX, the second faintest known dSph, 

has (v 2 ) 1/2 = 6.8±l;g (Chapman et al. 2005). 

Is there a minimum halo size in which stars can form, 
or a minimum clustering scale for the dark matter? There 
was no survivability reason why UMajor could not have 
been much less massive. For instance, the mass required 
to bind a 500 pc dwarf halo against tidal disruption by a 
10 12 Mq Galaxy at UMajor's distance of 100 kpc is es- 
timated by equating the dwarf and Galactic densities, so 
that Mbmd ~ 10 12 M & x (500pc/100kpc) 3 - 10 5 M Q , or 
about two orders of magnitude less massive than UMajor's 
central region, and perhaps three orders of magnitude less 
massive than its global mass. 

Kormendy & Freeman (2004) have compiled scaling re- 
lations among dwarf galaxies and more massive ellipticals, 
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relating Mb, po, r corc , and (v 2 ) ' . We can locate UMa- 
jor in this ensemble, using the values po = 0.18 M Q pc~ 3 , 
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(u 2 ) = 9.3 kms -1 , r coro <~ r 1 / 2 ~ 250pc, and Mb = 
My + 0.9 = —5.65, where My — —6.75 has been converted 
to B band using the typical B — V « 0.9 color of a K-giant. 
From Figures 2 and 4 of Kormendy & Freeman (2004), we 
find, unsurprisingly, that UMajor is typical of the dSphs in 
all respects but luminosity. All but two of the previously 
known dSphs and dlrrs fall onto a tight Faber- Jackson re- 
lation between velocity dispersion and luminosity, albeit 
with a systematic offset relative to larger ellipticals. UMi 
and Draco are the exceptions, being about three magni- 
tudes under-luminous for their mass. To agree with the 
general trend, UMajor would need to have a dispersion of 
< 3 kms -1 , or would have to be about five magnitudes 
brighter. 

It has been suggested that the dSphs' extreme M/L 
results from the expulsion of baryons from their shallow 
potential well by supernovae (Larson 1974; Dckcl & Silk 
1986), Because the potential wells of the dSphs seem com- 
parable, UMajor must have been intrinsically baryon poor 
to begin with, or else the baryon expulsion efficiency must 
have varied by one or two orders of magnitude among 
dSphs in order to produce their very different luminosi- 
ties. UMajor demonstrates that the luminosity scatter at 
the low mass end of the galaxy distribution is very large. 
It appears likely that more dark and massive dwarfs are 
lurking in the vicinity of the Galaxy. Detections of more 
candidate objects are urgently needed to check whether the 
number and properties of the population are truly consis- 
tent with the missing dark matter haloes of the simula- 
tions. 
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Fig. 1. — SDSS color-magnitude diagram for stars in the inner 6' of Ursa Major, with clearly visible giant and horizontal branches. The 
filled circles represent giant branch candidates observed in this work. 
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Fig. 2. — Velocities of Ursa Major candidate objects and <T SC at uncertainties. Stars 3 and 4 (near top) appear to be non-members, with 
velocities displaced about 60 kms -1 from the others. 
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Fig. 3. — Equivalent widths of the three Calcium triplet lines. The reddest two lines of objects 3 and 4, which are apparent non-members 
on the basis of velocity, also have a significantly larger equivalent width than the other stars. The equivalent widths for object 7 are more 
uncertain than the error bars indicate because the absorption lines were almost invisible to the eye, and the fits may be entirely spurious. No 
fit could be obtained for the reddest line of object 7. 



8 



Table 1 

Positions, magnitudes, and velocities of UMajor candidate stars 



Object 


RA 


Dec 


9 


i 


V 


0IRAF 


^scat 




(J2000) 






( kms -1 ) 


( kms -1 ) 


( kms -1 ) 


1 


10 h 35 m 28?51 


+51 57W9 


18.64 


17.45 


-54.87 


1.00 


1.30 


2 


10 h 34 m 30?51 


+51°57'07!0 


18.88 


17.75 


-50.80 


0.96 


1.09 


3 


10 h 35 m 15?87 


+51°59'32!0 


18.93 


17.82 


15.12 


1.45 


0.94 


4 


10 h 34 m 52M4 


+51°57'02{2 


19.08 


17.89 


8.21 


1.50 


0.40 


5 


10 h 34 m 52?05 


+51°58'28:S 


19.29 


18.25 


-64.22 


1.90 


0.65 


6 


10 h 34 m 42?36 


+51°58'06(1 


19.45 


18.37 


-38.95 


1.80 


1.25 


7 


10 h 35 m 17?23 


+51°55'33"7 


19.48 


18.38 


-53.53 


5.00 


5.15 



Note. - First column is an object identification number. Columns 2 and 3 are the position in J2000 coordinates. 
Columns 4 and 5 are the SDSS g and i magnitudes. Columns 6, 7, and 8 are v, the measured heliocentric velocity; ctiraf, 
the uncertainty obtained from the rescaled IRAF FXCOR velocity error; and cr sca t, the uncertainty obtained from the 
velocity scatter among the different absorption lines. 



